Correctness by construction for probabilistic programs

The "correct by construction" paradigm is an important component of modern Formal Methods, and here we use the probabilistic Guarded-Command Language $\mathit{pGCL}$ to illustrate its application to $\mathit{probabilistic}$ programming. $\mathit{pGCL}$ extends Dijkstra's guarded-command language $\mathit{GCL}$ with probabilistic choice, and is equipped with a correctness-preserving refinement relation $(\sqsubseteq)$ that enables compact, abstract specifications of probabilistic properties to be transformed gradually to concrete, executable code by applying mathematical insights in a systematic and layered way. Characteristically for "correctness by construction", as far as possible the reasoning in each refinement-step layer does not depend on earlier layers, and does not affect later ones. We demonstrate the technique by deriving a fair-coin implementation of any given discrete probability distribution. In the special case of simulating a fair die, our correct-by-construction algorithm turns out to be "within spitting distance" of Knuth and Yao's optimal solution

Value Iteration for Long-run Average Reward in Markov Decision Processes

Markov decision processes (MDPs) are standard models for probabilistic systems with non-deterministic behaviours. Long-run average rewards provide a mathematically elegant formalism for expressing long term performance. Value iteration (VI) is one of the simplest and most efficient algorithmic approaches to MDPs with other properties, such as reachability objectives. Unfortunately, a naive extension of VI does not work for MDPs with long-run average rewards, as there is no known stopping criterion. In this work our contributions are threefold. (1) We refute a conjecture related to stopping criteria for MDPs with long-run average rewards. (2) We present two practical algorithms for MDPs with long-run average rewards based on VI. First, we show that a combination of applying VI locally for each maximal end-component (MEC) and VI for reachability objectives can provide approximation guarantees. Second, extending the above approach with a simulation-guided on-demand variant of VI, we present an anytime algorithm that is able to deal with very large models. (3) Finally, we present experimental results showing that our methods significantly outperform the standard approaches on several benchmarks

Species-level functional profiling of metagenomes and metatranscriptomes.

Functional profiles of microbial communities are typically generated using comprehensive metagenomic or metatranscriptomic sequence read searches, which are time-consuming, prone to spurious mapping, and often limited to community-level quantification. We developed HUMAnN2, a tiered search strategy that enables fast, accurate, and species-resolved functional profiling of host-associated and environmental communities. HUMAnN2 identifies a community's known species, aligns reads to their pangenomes, performs translated search on unclassified reads, and finally quantifies gene families and pathways. Relative to pure translated search, HUMAnN2 is faster and produces more accurate gene family profiles. We applied HUMAnN2 to study clinal variation in marine metabolism, ecological contribution patterns among human microbiome pathways, variation in species' genomic versus transcriptional contributions, and strain profiling. Further, we introduce 'contributional diversity' to explain patterns of ecological assembly across different microbial community types

Relating circulating thyroid hormone concentrations to serum interleukins-6 and -10 in association with non-thyroidal illnesses including chronic renal insufficiency

<p>Abstract</p> <p>Background</p> <p>Because of the possible role of cytokines including interleukins (IL) in systemic non-thyroidal illnesses' (NTI) pathogenesis and consequently the frequently associated alterations in thyroid hormone (TH) concentrations constituting the euthyroid sick syndrome (ESS), we aimed in this research to elucidate the possible relation between IL-6 & IL-10 and any documented ESS in a cohort of patients with NTI.</p> <p>Methods</p> <p>Sixty patients and twenty healthy volunteers were recruited. The patients were subdivided into three subgroups depending on their underlying NTI and included 20 patients with chronic renal insufficiency (CRI), congestive heart failure (CHF), and ICU patients with myocardial infarction (MI). Determination of the circulating serum levels of IL-6 and IL-10, thyroid stimulating hormone (TSH), as well as total T4 and T3 was carried out.</p> <p>Results</p> <p>In the whole group of patients, we detected a significantly lower T3 and T4 levels compared to control subjects (0.938 ± 0.477 vs 1.345 ± 0.44 nmol/L, p = 0.001 and 47.9 ± 28.41 vs 108 ± 19.49 nmol/L, p < 0.0001 respectively) while the TSH level was normal (1.08+0.518 μIU/L). Further, IL-6 was substantially higher above controls' levels (105.18 ± 72.01 vs 3.35 ± 1.18 ng/L, p < 0.00001) and correlated negatively with both T3 and T4 (r = -0.620, p < 0.0001 & -0.267, p < 0.001, respectively). Similarly was IL-10 level (74.13 ± 52.99 vs 2.64 ± 0.92 ng/ml, p < 0.00001) that correlated negatively with T3 (r = -0.512, p < 0.0001) but not T4. Interestingly, both interleukins correlated positively (r = 0.770, p = <0.001). Moreover, IL-6 (R²= 0.338, p = 0.001) and not IL-10 was a predictor of low T3 levels with only a borderline significance for T4 (R²= 0.082, p = 0.071).</p> <p>By subgroup analysis, the proportion of patients with subnormal T3, T4, and TSH levels was highest in the MI patients (70%, 70%, and 72%, respectively) who displayed the greatest IL-6 and IL-10 concentrations (192.5 ± 45.1 ng/L & 122.95 ± 46.1 ng/L, respectively) compared with CHF (82.95 ± 28.9 ng/L & 69.05 ± 44.0 ng/L, respectively) and CRI patients (40.05 ± 28.9 ng/L & 30.4 ± 10.6 ng/L, respectively). Surprisingly, CRI patients showed the least disturbance in IL-6 and IL-10 despite the lower levels of T3, T4, and TSH in a higher proportion of them compared to CHF patients (40%, 45%, & 26% vs 35%, 25%, & 18%, respectively).</p> <p>Conclusion</p> <p>the high prevalence of ESS we detected in NTI including CRI may be linked to IL-6 and IL-10 alterations. Further, perturbation of IL-6 and not IL-10 might be involved in ESS pathogenesis although it is not the only key player as suggested by our findings in CRI.</p

Analysing RoboChart with probabilities

Robotic systems have applications in many real-life scenarios, ranging from household cleaning to critical operations. RoboChart is a graphical language for describing robotic controllers designed specifically for autonomous and mobile robots, providing architectural constructs to identify the requirements for a robotic platform. It also provides a formal semantics in CSP. RoboChart has a probabilistic operator (P) but no associated probabilistic CSP semantics. When (P) is used, currently a non-deterministic choice (Π) is used as semantics; this is a conservative semantics but it does not allow the analysis of stochastic properties. In this paper we define the semantics of the operator in terms of the probabilistic CSP operator ⊞. We also show how this augmented CSP semantics for RoboChart can be translated into the PRISM probabilistic language to be able to check stochastic properties

MiR-100 regulates cell differentiation and survival by targeting RBSP3, a phosphatase-like tumor suppressor in acute myeloid leukemia

Acute myeloblastic leukemia (AML) is characterized by the accumulation of abnormal myeloblasts (mainly granulocyte or monocyte precursors) in the bone marrow and blood. Though great progress has been made for improvement in clinical treatment during the past decades, only minority with AML achieve long-term survival. Therefore, further understanding mechanisms of leukemogenesis and exploring novel therapeutic strategies are still crucial for improving disease outcome. MicroRNA-100 (miR-100), a small non-coding RNA molecule, has been reported as a frequent event aberrantly expressed in patients with AML; however, the molecular basis for this phenotype and the statuses of its downstream targets have not yet been elucidated. In the present study, we found that the expression level of miR-100 in vivo was related to the stage of the maturation block underlying the subtypes of myeloid leukemia. In vitro experiments further demonstrated that miR-100 was required to promote the cell proliferation of promyelocytic blasts and arrest them differentiated to granulocyte/monocyte lineages. Significantly, we identified RBSP3, a phosphatase-like tumor suppressor, as a bona fide target of miR-100 and validated that RBSP3 was involved in cell differentiation and survival in AML. Moreover, we revealed a new pathway that miR-100 regulates G1/S transition and S-phase entry and blocks the terminal differentiation by targeting RBSP3, which partly in turn modulates the cell cycle effectors pRB/E2F1 in AML. These events promoted cell proliferation and blocked granulocyte/monocyte differentiation. Our data highlight an important role of miR-100 in the molecular etiology of AML, and implicate the potential application of miR-100 in cancer therapy

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO's second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h095%=3.47×10-25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering. © 2019 American Physical Society

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO’s second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h95%0=3.47×10−25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering